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Abstract. The volume coherent deflection of high-energy positive and negative 
particles in uniformly bent crystals is studied. The general analysis of potential 
scattering shows that the standard screening potential for a large class of collisions 
can cause the volume refraction for negative particles (antiprotons, electrons) 
instead of the volume reflection for positive particles (proton, positrons). 



PACS numbers: 61.85.+P; 45.50Tn; 41.85-p; 25.75.Gz 

1. Introduction 

In the early 80-s, the experimental studies [11 [21 [3] have discovered the volume capture 
of relativistic particles in a channeling motion by bent crystals. Several attempts 
have been made to find a specific mechanism of capture(see e.g.[i4j) because the elastic 
scattering in a central field does not include the channeling trajectories (see below). 
The computer simulations of the passage of high energy particles in a uniformly bent 
crystal [3 [6] were completed and it was found that the particles (both positive and 
negative) are deflected to the side opposite the bend of atomic planes by an angle of 
about 2 * 6'c {9c - critical angle). It was called a "volume reflection" by bent atomic 
planes. The signiflcant beam divergence and multiple scattering of the quasi-channeled 
particles in the crystals did not allow to find this effect for almost two decades. The 
experiment [7] [8] of lOGeV protons scattered by a short bent crystal is probably the 
first explicit demonstration of the volume reflection. Early the preliminary results of 
crystal collimation at RHIC J, for the 2001 run demonstrated some interested features 
which were later explained [TUl [TT] by the same effect. 

Our purpose is to discuss the principle difference between deflections of positive 
and negative particles if the standard continuous screening potential for the crystal 
planes is considered to be a good approximation. In this approximation the sign of 
continuous screening potential is deflned by the charge of projectile and the behavior 
of positive and negative particles are different. The positive particles should undergo 
the "volume reflection", but it turns out that the negative particles in such potential 
should deflect in the same direction as the bend of atomic planes. It is not the 
channeling motion because the elastic scattering in central potential produces two 
symmetrical branches touching equipotential circle and extending from the turning 
point to infinity. This effect is different from predicted in [51 [S] and can be called a 
"volume refraction" by the bent atomic planes. 
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2. Screening Crystal Potential 

The charged particles incident on a single crystal with small angles to the 
crystallographic directions, experience the collective fields from atoms of slightly bent 
planes. The interaction potential for binary ion-atom collisions is, in general, described 

by 
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where Zp^T are the bare nuclear charge of the projectile and crystal atoms, Rt-f is 
Thomas-Fermi (or Firsov) screening length and / is a screening function (/ > 0). The 
screening function is 
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and it approaches 1 in the limit r 
only possible if 



Rt-f 

and decays exponentially for large r. This is 



= 1. 



Different analytic forms of the screening function / are in use. The Moliere 
approximationUg, /a/, uses {a,} = {0.1,0.55,0.35}, {/?,} = {6.0,1.2,0.3}. Other 
choices include the Ziegler-Biersack-Littmark potential[13 fzBL with {at} = 
{0.1818, 0.5099, 0.2802, 0.0281}, {/3^} = {3.2,0.9423,0.4028,0.2016}. Now, if we find 
the average potential over a cylinder with the radius R » Rt-f and the root-mean- 
square temperature displacement u of atoms in gaussian form, we receive the following 
expression for continuous cylindrical potential 
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Here p is polar coordinate from center of cylinder, U — ttZpZtRt-fit-p, Up is a density 
of atoms in the bent plane, erfcQ is the complementary error function|14j. 

U 



Figure 1. The radial profile of cylindrical potential. 

For positive {Zp > 0) or negative {Zp < 0) charged projectiles , the screen 
potential of atoms H]) over all space (0 < r < oo) is repulsive {Uatomir) > 0) or 
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attractive {Uatom{f) < 0) respectively. (Strictly speaking, the potential is called 
repulsive if the force is everywhere repulsive, i.e., the potential monotonely decreases 
to zero, but for simplicity we use the above terminology.) It follows that the average 
potential of each plane Up{p) and the whole cylindrical crystal potential 



will also be repulsive or attractive depending on the charge of the projectile (Fig. 
[ij. The summation in ^ is over a range of radii Ri,i = 1..N which denote the 
centers of atomic planes located periodically from Ri to Rn- For < p « Ri and 
Rn « p < oo the potential decreases to zero. In the area Ri < p < Rm, potential 
does not have a singularity and is smooth almost periodic function with some constant 
positive or negative pedestal. The value of pedestal is quite significant. It is 13.3 
eV for Si < 110 > plane (T=300K), i.e. about 1/3 of of maximal value of screening 
potential ( 35.9 eV in this case). 

3. Scattering in Central Field 

The path of a particle in a central classical field is, as well known, symmetrical about a 
line from the center to the nearest point in the orbit (OA for positive and OB negative 
particles in Fig[2]). Hence the two asymptotes to the orbit make equal angles {(p+ for 



positive and for negative particles) with these lines. The angles x+.- = ti" ~ "2^+.- 
through which the positive and negative particles are deflected as they pass the center 
of potential have the different sign (see Fig[21), which we discuss below. This general 
rule is applied to any centrally (cylindrically) symmetrical potential including the 
cylindrical crystal with radius R. 

The solution for relativistic Hamilton-Jacobi equation in central field can be found 
( see e.g. [15]) and it gives the well known expression for tp: 




(4) 




Figure 2. Scattering by repulsive and attractive potentials 
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where E = \/'p\^(? + m^c**, lJ{p) are the energy and potential of central field 
^, Pmim Poc are the impact parameter, the nearest point in the orbit (p„ji„=OA or 
Pmin—OJi in Figl2|) and the momentum of particle at infinity respectively. If we use 
new variable g = p/b, the scattering angle ^ can be written in more convenient form: 

where g{g) = '^'^■i'^l^ . We should note that similarity of g{g) and the square of 
Lindhard angle 0^ — is not accidental. The nearest point Qmin is defined by the 

Poo ^ 

transcendental equation 



(7) 
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In the case of purely attractive potential, g{g) < 0, from ([7]) we can receive gmin < 1, 
i.e. the nearest point gmin is always less than impact parameter. In the case of purely 
repulsive potential g{g) > 0, we have opposite relation gmin > 1- In the case of 
U (bg) = 0, the nearest point gmin = 1 and the integral ^ becomes 

which means that the particle does not deflect from straight line. 
4. Classical Limits of Scattering Angle 

Now we can find the limits for classical variation of scattering angles ip^ , y+ . As it 
was pointed out for attractive potential < 0), Pmin is less than impact parameter 
b and the integral ([6]) can be estimated as 
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Here g_ = max{\g{g)\), gmin- < < oo is a deepest value of the attractive potential 
(H]). The (HI) gives the lower bound of the scattering angle (p^. For positive particles 
> 0), it is easy to see that 
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where g+ = max{g{g)) > g{g) is the highest point of repulsive potential. This is 
different from the scattering on repulsive monotonely decreasing potential [TS] , where 
this angle is always less than 7r/2. 
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5. Small Angle Scattering 

Now let us turn to the small angle scattering in the central field ([4]). From physical 
pictures of scattering (Figl2]), it is clear that the particle's path has a biggest curvature 
near the turning point. So, the potential near that point plays a crucial role in the 
forming of scattering angles. In the rest of space the potential is analytic function, 
so we can introduce 'residual' potential g — g{Qmin) — 9- Since g is small in the 
neighborhood of Qmim we expand the square root \J {1 — 9{Qmin) + 'g)Q^ — 1 (eq. (O) 
in powers of this residual potential: 
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The first integral gives 7r/2. The second integral depends on the magnitude and 
sign of scattering 'residual' potential, as well as on the impact parameter and energy 
of particles. Its integrand has a divergent term in denominator at g — Qmin- However, 
the 'residual' potential 'g{gmin) — 0, and there is no divergence in integral. The second 
integral can be integrated by parts, giving 

= 1-- n ^ nn 

Finally, using || = — ^ the formula for classical deflection angle for particles is 
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The difference between deflection of positive and negative particles is seen from the 
following equations: 
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6. Conclusion 



The angle of scattering given by these equations is the first approximation in the 
expansion pip in powers of the residual potential and can describe the major features 
of small angle deflection of positive and negative particles. It is more general than the 
classical small angle formula (see e.g. [16]) in sense that it describes more accurate the 
turning points. They coincide if potentials have a long range of transversal variations 
and £imm_,+ — > 1, 5(f?mm_ +) << 1 in (|13|14|) . In this case we have 

f°° dg dg 

''^^-IjgW^' ^''^ 
, ^9 dg 



A POSSIBILITY OF VOLUME REFRACTION OF NEGATIVE RELATIVISTIC PARTICLES IN BENT CRYSTALS 



Since U+{p) — —U^{p), the deflective angles for positive and negative particles 
satisfy the following conditions: 



X+ = -X- (17) 

So if the positive particles deflected by positive angle x+, the negative particles 
should deflected to the opposite direction. For more accurate calculations one should 
use the formulas (|13ll4p . 
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